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A new scheme of field quantization is proposed. Instead of associating with different frequencies 
different oscillators we begin with a single oscillator that can exist in a (quantum) superposition of 
different frequencies. The idea is applied to the electromagnetic radiation field and nonrelativistic 
quantum optics. Employing a Dirac-type mode-quantization of the electromagnetic field and using 
a single oscillator we obtain several standard properties such as coherent states or spontaneous 
and stimulated emission. Extending the formalism to a greater number of oscillators we arrive 
at a structure analogous to the Fock space but without the standard cyclic "vacuum state". In 
the modified formalism the notion of the vacuum state is replaced by a vacuum subspace spanned 
by ground states of the oscillators. As opposed to the standard approach the vacuum energy is 
finite and does not have to be removed by any ad hoc procedure. Atom-light interactions are 
described by an appropriately modified minimal-coupling Hamiltonian (no normal ordering of the 
free-field Hamiltonian is necessary). The Hamiltonian does not change the number of oscillators 
which leads to an additional conservation law. Using the "— (e/m)A ■ p" interaction we discuss in 
second-order perturbation theory a two-photon spontaneous emission. The result essentially agrees 
with the ordinary formulas but the nontrivial vacuum structure is explicitly seen in the two-photon 
amplitude. The probability of the 2-photon emission resulting from the new formalism consists of 
a product of several terms, a part of them resembling those arising in the standard formulation 
from detector inefficiency, and the remaining one being the well known quantum optics formula. 
The presence of the additional conservation law shows that the theory contains two kinds of bosons 
(oscillators, whose number is conserved, and their excitations, whose number is not conserved). 
Taking this distinction into account we calculate an analog of the blackbody radiation Planck law. 
For temperatures lower than some T cr iticai the result is indistinguishable from the Planck distribution. 
For T > T C ri t i ca i the distribution is Planck-like but with the maximum lowered and shifted towards 
higher frequencies. 



I. INTRODUCTION 



The standard quantization of a harmonic oscillator is based on quantization of p and q but w is a parameter. To 
have, say, two different frequencies one has to consider two independent oscillators. On the other hand, it is evident 
that there can exist oscillators which are in a quantum superposition of different frequencies. The example is an 
oscillator wave packet associated with distribution of ccntcr-of-mass momenta. It is known that the superposition of 
momenta gets translated into a superposition od Dopplcr shifts and therefore also of frequencies. We stress here the 
word "quantum" since the superpositions we have in mind are not those we know from classical oscillations. 

This trivial observation raises the question of the role of superpositions of frequencies for a description of a single 
harmonic oscillator. The motivation behind the problem is associated with the question of field quantization: Is it 
possible that a quantum field consists of oscillators whose frequencies are indefinite? If so, maybe to quantize the 
field it is sufficient to use only one oscillator which exists in a quantum superposition of all the possible frequencies 
allowed by the boundary conditions of a given problem? 

The idea is very simple. It is known that a "one-particle" state vector can be regarded as a representation of an 
ensemble of particles in a given pure state. On the other hand, the classical electromagnetic field can be regarded 
as an ensemble of oscillators. The standard idea of quantization, going back to 1925 W, is to treat the field as an 
ensemble of quantum oscillators. But the ensemble itself is, in a sense, a classical one since for each frequency we need 
a separate oscillator. This is analogous to a classical ensemble of particles forming a classical wave on a lake surface. 
For each point on the surface we need a separate particle because a classical particle can ocupy only a single point in 
space. A quantum wave is of course different and we are all accustomed to the idea of a single-particle wave. In this 
case the properties of the entire ensemble are somehow encoded in properties of a single element of the ensemble. 

For some reasons, probably partly historical and sociological, it seems that the idea of a single-particle state vector 
representation of the ensemble of oscillators has never been considered. The historical reason may be the fact that 
the very concept of field quantization occured already in 1925. At that stage quantum mechanics existed still in a 



1 



matrix form and the Schrodinger paper "Quantisierung als Eigenwertproblem" |2j, where the Schrodinger equation 
occured for the first time and the role of eigenvalues was explained, was not yet published. Sociologically, the names 
and reputation of Heisenberg, Born, Jordan, Dirac, together with the unquestionable success of quantum optics, field 
theory, and statistical physics, made it almost impossible to question the very starting point of the theory. The ideas 
presented below are an accidental by-product of a work on a different problem. 

It should be mentioned that several approaches towards an alternative description of the electromagnetic field at 
a fundamental level were proposed (e.g. Janes' || neoclassical theory, stochastic electrodynamics [||). But the main 
idea of all such alternatives was to treat the field in classical terms and to associate the observed discreteness of 
emission/absorbtion phenomena with the quantum nature of atoms and not with the field itself. 

The approach we will discuss in this paper does not belong to this tradition, is much more radical and goes in the 
opposite direction. We will not try to make the field more classical. What we will try to do is to make it even more 
quantum by replacing classical parameters with eigenvalues. 

The field will be quantized at a one-particle level, but then extended to multi-particle systems. Obviously, it is not 
possible to include in a single paper all the possible tests of the new formalism one should perform. We will therefore 
concentrate on these points where quantum electrodynamics produces results which are believed to be a consequence 
of the standard canonical quantization. Three areas should be checked first: 

(i) Vacuum effects in atomic physics. 

(ii) Emission of photons in entangled states. 

(iii) Boson statistics and the Planck law. 

For the first two problems we shall choose the simplest approach, namely first and second order perturbation theory. 
The Planck law will be discussed in a more detailed way. As we shall see the new theory is not completely equivalent 
to the standard one, but the modifications one finds are surprisingly subtle and in principle subject to experimental 
tests. 

In a separate paper we shall discuss perturbation theory to all orders, since then a kind of prescription for translating 
the old results into the new framework may appear. 



II. HARMONIC OSCILLATOR IN SUPERPOSITION OF FREQUENCIES 

We know that frequency is typically associated with an eigenvalue of some Hamiltonian or, which is basically the 
same, with boundary conditions. A natural way of incorporating different frequencies into a single harmonic oscillator 
is by means of the frequency operator 

0= ^2 u k \uJk,jk)(uJk,jk\ (1) 

where all ujk > 0. For simplicity we have limited the discussion to the discrete spectrum but it is useful to include 
from the outset the possibility of degeneracies, represented here by the additional discrete quantum numbers jk- The 
corresponding Hamiltonian is defined by 

H = Ml ®-(a t a + aa f ) (2) 
where a = X^J^Lo V n + l\n)(n + 1|. The eigenstates of H are \u)k,jk,ri) and satisfy the required formula 

H\uj k ,j k ,n) = huj k (n+ ^j\aj k ,jk,n) (3) 

justifying our choice of H. The standard case of the oscillator whose frequency is just ui coresponds either to = loI 
or to the subspace spanned by \LOk,jk,n) with fixed u>k = lo. Introducing the operators 

a u k ,jk = l w fe>Jfe)< w fe)Jfel ® a (4) 

we find that 

H =\Yj n ^( a i k ,3 k a ^» + a ^^ a L, ]k )- (5) 
^>k -Jk 

The algebra of the oscillator is "noncanonical" : 
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Kfcj^a^jJ = ^ k ^S jkjl \uj k ,j k )(uj k ,j k \ ® 1 (6) 

a u k ,j k a u)i,ji — S L j kU)l 6j k j l (a UJk j k ) (7) 

a L,jk a liJ, = 6 "k"l S Ml( a u k ,j k f- ( 8 ) 

The dynamics in the Schrodinger picture is given by 

tBft|*> = fT|*) = hSl 0(0*0 + (9) 
In the Heisenberg picture we obtain the important formula (see Appendix XI A) 

a. kljk (t)=e iHt/R a Ukdk e- iHt / h (10) 

2 ^k ,jk ■ 

Taking a general state 



|w fc ,ifc)(wfc,ifc| 8e-" kt a = e^'o^,. (11) 



X] ^(^fe, Jfc,«)l^fc, jfc)l"-> ( 12 ) 



we find that the average energy of the oscillator is 

(H) = (ip\H\ip) = £ |#^,j fe ,n)| 2 7ia; fe (n+ J). (13) 



The average clearly looks as an average energy of an ensemble of different and independent oscillators . The ground 
state of the ensemble, i.e. the one with tjj(uj k ,j k , n > 0) = has energy 

(H) =\ ]T |^K,j fc ,0)| 2 ^ fe (14) 

U> k ,jk 

which is finite if 

^2 "^(^kJkiO^uJkJk) (15) 

belongs to the domain of $7. The result is not surprising but still quite remarkable if one thinks of the problem of 
field quantization. 

The very idea of quantizing the electromagnetic field, as put forward by Born, Heisenberg, Jordan []l] and Dirac ||, 
is based on the observation that the mode decomposition of the electromagnetic energy is analogous to the energy of 
an ensemble of independent harmonic oscillators. In 1925, after the work of Heisenberg, it was clear what to do: One 
had to replace each classical oscillator by a quantum one. But since each oscillator had a definite frequency, to have 
an infinite number of different frequencies one needed an infinite number of oscillators. The price one payed for this 
assumption was the infinite energy of the electromagnetic vacuum. 

The infinity is regarded as an "easy" one since one can get rid of it by redefining the Hamiltonian and removing 
the infinite term. The result looks correct and many properties typical of a quantum harmonic oscillator are indeed 
observed in electromagnetic field. However, subtraction of infinite terms is in mathematics as forbidden as division by 
zero so to avoid evident absurdities one is forced to invent various ad hoc regularizations whose only justification is 
that otherwise the theory would not work. In larger perspective (say, in cosmology) it is not at all clear that an infinite 
(or arbitrarily cut off at the Planck scale) energy of the vacuum does not lead to contradictions with observational 
data ||. Finally, Dirac himself had never been fully satisfied by the theory he created. As Weinberg put it, Dirac's 
"demand for a completely finite theory is similar to a host of other aesthetic judgements that theoretical physicists 
always need to make" . 

The oscillator that can exist in superpositions of different frequencies is a natural candidate as a starting point for 
Dirac-type field quantization. Symbolically, if the Heisenberg quantization is p 2 + oj 2 q 2 i— > p 2 + uj 2 q 2 , where w is a 
parameter, the new scheme is p 2 +uj 2 q 2 i— ► p 2 -\-Cj 2 q 2 , where to is an operator. Its spectrum can be related to boundary 
conditions imposed on the fields. The field now can exist in superposition of frequencies but the superposition is meant 
in the quantum sense i.e. the field may consist of (an indefinite number) of oscillators with indefinite frequency. In 
this meaning the approach we propose is even "more quantum" than the standard one since u> is not a (classical) 
parameter but an eigenvalue. 

We do not need to remove the ground state energy since in the Hilbert space of physical states the correction is 
finite. The question we have to understand is whether one can obtain the well known quantum properties of the 
radiation field by this type of quantization. 
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III. PRELUDE: "FIRST QUANTIZATION" — FIELD OPERATORS FOR FREE MAXWELL FIELDS 

The new quantization will be performed in two steps. In this section we describe the first step, a kind of first 
quantization. In next sections we shall perform an analogue of second quantization which will lead to the final 
framework. It is essential that the "second quantization" will not involve, in fact, any additional quantization but is 
simply a transition from one to many oscillators. 

The energy and momentum operators of the field are defined in analogy to H from the previous section 

H = ^2 hu\\s,K\){s,K\\ ® ^{a^a + aa^ (16) 

S,K A 

= \ Yl hUJ *( a t,Kx a ^X + fl S,«AClL A ) (17) 
S,K A 

P= ^Mx\s,K\){s,Kx\®^{a { a + aa [ ^ (18) 

S,K A 

= \ Yl ^("La^KA + a S ,KA a L A ) ( 19 ) 



S,K A 



where s = ±1 corresponds to circular polarizations. Denote P — {H/c, P) and P ■ x = Ht — P ■ x. We employ the 
standard Dirac-type definitions for mode quantization in volume V 



A(t, f) = (a s ^e-^e s ^- s + a^^'e^e"^) (20) 

= e iP^/^Ae- lP - x/n (21) 
^) -^^{W^^^ < KA ^V^e* KA ) (22) 

S,K A 

= e iP-*/KE e -iP-x/t> (23) 

(24) 



(25) 



where 



a s ,« A = \s,kx)(s,k\\ ®a (27) 
al z x = \s,k x )(s,kx\ <W- (28) 

For later purposes we introduce the notation 

[a s ,K A ,a^ s J = 1 s ,k a = \s,K\)(s,FZ\\ <g> 1. (29) 
Now take a state (say, in the Heisenberg picture) 

I*) = E y s ,K X ,n\s,K X ,n} (30) 

= E $ s,kJ S >«a)K,k a ) (31) 

S,K A 

where |a s ,re A ) form a family of one-oscillator coherent states: 

a\a s ,z x ) = a s ^ x \a Si u x ) (32) 
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The averages of the field operators are 

(*|l(t,f)|*) = ]T|$ s ^ C!:!! 



*|i(t,f)|*> = J2 !*-,»> I V^^x^-"^^ -< KA e->X«J ( 34 ) 



=i^|$^ A |y^(a s>K ,e-^-n K , x e s , KA - < KA e lK ^n KA x e s * K J (35) 

These are just the classical fields. More precisely, the fields look like averages of monochromatic coherent states with 
probabilities I^kJ 2 . The energy-momentum operators satisfy also the standard relations (see Appendix XI B) 



H =\J d 3 x^E(t,x)-E(t,x)+B(t,x)-B(t,x)y (36) 
P= I d 3 xE(t,x) x B(t,x). (37) 



Jv 

It should be stressed, however, that these relations have a completely different mathematical origin than in the 
usual formalism where the integrals are necessary in order to make plane waves into an orthonormal basis. Here 
orthogonality follows from the presence of the projectors in the definition of a SiKx and the integration in itself is trivial 
since 

E{t, x) ■ E{t, x) + B(t, x) ■ B(t, x) — E ■ E + B ■ B (38) 
E(t, x) x B(t, x) = E x B. (39) 

Therefore the role of the integral is simply to produce the factor V which cancels with 1/V arising from the term 
1/VV occuring in the mode decomposition of the fields. To end this section let us note that 

(*|ff|*> = ^ A |$ S , KA | 2 (K KA | 2 + i) (40) 

= J2 ^a|$ s , ka i 2 (Wj 2 + 1). (41) 

S,K X 

The contribution from the vacuum fluctuations is nonzero but finite. One can phrase the latter property also as 
follows. The noncanonical algebra of creation- annihilation operators satisfies the resolution of identity 

EKA.«U] = 1 (42) 

wheras the canonical algebra would impliy 

EK*»' a U] =ooL (43) 

IV. SPONTANEOUS AND STIMULATED EMISSION: FIRST VERSION 

Some typically quantum optical phenomena occur already at the one-oscillator level. Below we shall see that 
spontaneous and stimulated emissions are a property of a single-oscillator description, although to have a more 
complete picture we need the multi-oscillator extension discussed in subsequent sections. 

Beginning with the dipole and rotating wave approximations (RWA) we arrive at the Hamiltonian 

H = ^huj a 3 + i huJx ( a l,n x a sM>, + a s ^A a !,K A ) + ^od ^ (gs,K X a s M x <7+ + 0* ^fl^ou) (44) 
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where du = (+\d\— ) is the matrix element of the dipole moment evaluated between the excited and ground states, 



and g s M x = i y 2hL> x v ^ s -Z* ' ^ ^ ne Hamiltonian represents a two-level atom located at xo = 0. 
The Hamiltonian in the interaction picture has the well known form 

Hj = huJodY, (g s ^"-^ t a s ^a + + g^e^^K^a-). (45) 
Consider the initial state 

= *'',%,aW>%>A+)+ * S V;,nl s '>'C^+>- (46) 

S ,K q s , tz ^1 

The states corresponding to n = play a role of a vacuum. As a consequence the vacuum is not represented here by 
a unique vector, but rather by a subspace of the Hilbert space of states. Energy of the general vacuum state 

l*) = ]T *.,* A) o,±l*>*A>0.±> (47) 

S,K A ,± 

is related to the density of modes p(k\) = J2 S ± |* a ,u A ,o,± | 2 an d is, therefore, state dependent. 

In order to estimate the probabilities of spontaneous and stimulated emissions we can use the first-order time- 
dependent perturbation theory || and arrive at 

|*(f)) = |*(0)) 

e -i(a;o-wx )i _ I 

+UJ ° d 2^ 7 7 * s >«a ,o9t,z Xo \s, k Xo , 1, -) 

S,K 

E e -t(u -UA„)< _ I 
y a ,H n ,nVn + lg* s zjs, K n , n + 1, -). (48) 

One recognizes here the well known contributions from spontaneous and stimulated emissions. It should be stressed 
that although the final result looks familiar, the mathematical details behind the calculation are different from what 
we are accustomed to. For example, instead of 

a li,Hi I s ' ^' TO ) ~ 1; s, k, m), (49) 

which would hold in the standard formalism for k\ ^ k, we get simply 

al iA \s,K,m)=0, (50) 

a consequence of of - = 0. 

Let us now look more closely at spontaneous emission (we take n = 0). The state vector is 

= |*(0)) 

e -i(u -ux )t _ i 

+u; dV Gin |s,ka„,1,-> (51) 

S,K,Q 

where we have introduced the "effective coupling terms" 

G*s,^ =^,H Xo ,ogU Xo . (52) 

As we can see the result is mathematically equivalent to the standard one but with the coupling constants automatically 

regularized by the presence of the vacuum amplitude in G* s . 



We apply the standard argument but with g's replaced by G"s and obtain the probability of emission per time unit 



G 



P = 2i:u; 2 d 2 \^sM,o9s,k\ 2 S(uJo - ur)- 

S,K 



(53) 



Assuming for simplicity that density of vacuum modes is isotropic and polarization independent we can write it as a 
function of frequency only i.e. 

|* s ^ | 2 - F{ug) (54) 

and finally 

P = 2™ 2 d 2 F{w ) \9sM 2 S(^o - wa) = F(u )P o u- (55) 

Here P id is the emission rate obtained in the standard theory. The nontrivial structure of the vacuum influences the 
lifetime of the atom. We shall return to this and related questions later but first have to extend the formalism in a 
way allowing to consider entangled states of light. 

More reliable estimates of the lifetime require more detailed calculations that we postpone to a forthcoming paper. 
One should also keep in mind the possibility of an unisotropic vacuum caused by more complicated boundary conditions 
such as those occuring in measurements of the Casimir force. 



V. "SECOND QUANTIZATION" 

The Hilbert space of states of the field we have constructed is spanned by vectors \s,K,n). Still there is no doubt 
that both in reality (and the standard formalism) there exist multiparticle entangled states such as those spanned by 
tensor products of the form 

| + ,Kl,l)|-,K 2 ,l>, (56) 

and the similar. It seems that there is no reason to limit our discussion to a single Hilbert space of a single oscillator. 
What we have done so far was a quantization of the electromagnetic field at the level of a "one-particle" Hilbert space. 
Similarly to quantization of other physical systems the next step is to consider many (noninteracting) particles. 

The procedure is essentially clear. Having the one-particle energy-momentum operators P a (i.e. generators of 
4-translations in the 1-particle Hilbert space) we define in the standard way their extensions to the Fock-type space 

Va = Pa 

©(P ®1 + 1®P ) 

©(P a ® 1 ® 1 + 1 ® P a ® 1 + 1 ® 1 ® P Q ) 

© . . . . (57) 
The ^-dependence of fields is introduced similarly to the one-particle level 

f(t, x) = ^-P-x/h^ e -iv-x/n ( 58 ) 

but the field itself has yet to be defined. Assume 

T = Cl F 

©c 2 (P ® 1 + 1 ® F) 

©c 3 (F® 10 1 + 1® F® 1 + 1® 1®F) 

e . . . (59) 

where Ck are constants discussed below, and F is A, E, or B. The multi-oscillator annihilation operator associated 
with such fields must be therefore of the form 

®c 2 (a s ^ ® 1 + 1 ® a s ^) 

ec 3 (a 5iS ® 1 ® 1 + 1 ® a s ^ ® 1 + 1 ® 1 ® ci s m) 

© . . . . (60) 
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Having two 1-particle operators, say X and Y, one can easily establish a relation between the 1-particle commutator 
[X, Y] and the commutator of the extensions X, y-. 

®c 2 2 ([X, y] ®1 + 1® [X, Y]) 

®cl ([x, y] <g> i ® i + 1 <g> [x, y] <g> i + 1 ® i <g> [x, y]) 

©.... (61) 
The annihilation operators so defined satisfy therefore the algebra 

K«,4,«']=° for( S ,«)^( S ',«'), (62) 

[°s,a> = ( 63 ) 

[o a ,R,O a ^] = (64) 

[<g.<«-]=0 (65) 

where the operator l s ^ is defined by 

1 - - r 2 1 - 



3c|(l s>S (g) 1 + 1®1 8)S ) 

3c|(l S) £ (g) 1 (g) 1 + 1 (g) l s> g ®l + l®l(gl l a>B ) 



(66) 



and l S) g is a single-oscillator operator ( p9| ) 
As opposed to the single-oscillator case 



a s ,F-a s ,^ r ^ S ss <6z^(a S: z) 2 . (67) 



An important property of the 1-oscillator description was the resolution of identity (42). The requirement that the 
same be valid at the multi oscillator level leads to c„ = \l Jn. In such a case one finds that 



1* a * l.,a (68) 

but nevertheless 

ZX* = 1 - (6 9 ) 

Below we shall give another justification of this particular choice of c n . 
We can finally write 



X{t,x) = y^^^^'^e^^ + o^e^'^e-^- 2 ) (70) 

= jV-x/n^r-iV-x/h ( 71 ) 



= i £ ^^(a s>K ,e-^V*-*X^ " at^e^^e-^e^J (72) 

(74) 

= e* p "/ A $e-* p " /A . (76) 

These operators form a basis of the modified version of nonrelativistic quantum optics. 

Let us return for the moment to the case of a general c„. A straightforward calculation shows that 



8 



J V S,K X 

(Bel ({a S; £, aj -} ® 1 + 1 ® {a s ^, a^-} + 2a\ R ® a S; £ + 2a S; £ <g> 
©c^ ({a s ,g, aj -} ® 1 ® 1 + 2a S;S ® aj - ® 1 + 2a Si g 1® 
+ 2aJ - ® a s , s 01 + 18 {a s ,«, a^} 1 + 21 a StR ® a\ R 
+ 2a\- 1® a s ^ + 21 aj - a sfi + 10 10 {a s ^, a[ -}) 



Comparing this with the generator of time translations 



)({a S; £, aj -} 1 + 10 {a S;S , »^}) 

al,g} 10 1 + 10 {a a , K , a\ ~} 1 + 1 



K,a>aJ,a}) 



(77) 



(78) 



(79) 



we can see that there is a relation between H and H but the latter contains terms describing interactions between the 
oscillators. The contribution from these interactions vanishes on vacuum states. Below, when we introduce the notion 
of a generalized coherent state, we will be able to relate averages of Ti and H. In a similar way one can introduce the 
"Pointing operator" 



-I 



d 3 x£(t,x) x B(t,x) 



S,K X 



(80) 



Its relation to the generator of 3-translations V is similar to this between H and H. 

In the above construction the only clement which is beyond a simple transition to many oscillators is the choice of 
c„. For different choices of these constants we obtain different algebras of noncanonical commutation relations and 
therefore also different quantization schemes. Several different ways of reasoning lead to c n = l/^/n as we shall also 
see in the next sections. 



VI. SOME PARTICULAR STATES 

We assume that all the multi-oscillator states are symmetric with respect to permutations of the oscillators. 



A. Generalized coherent states 



For general c n an eigenstate of a s ,K A corresponding to the eigenvalue a St z is of the form 

|« s ,«) = fi(s, k)\s, k, a s ^/ci) 

®/ 2 (s, k)\s, k, a St nl (2c 2 )) \s, k, a St f>/ (2c 2 )) 

®h(s, «)k k, a s fi/ (3c 3 )) \s, k, a s fi/ (3c 3 ))|s, k, a s ^j (3c 3 )) 



(81) 
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where 

\s,n,a SiH ) = \s,n)\a SiH ), (82) 

J2k \fk(s,K.)\ 2 = 1) an( i a \ a s,R) — ots,n\o<-s,R)- What is interesting not all fk have to be nonvanishing. 
The average "energies" of the field in the above eigenstate are 



1 | 

(a sM \n\a s ,K) = ft^las.gp^-r-j |/ fc (s,K)| 2 + -Tilo r E k\f k (s, k)\ 2 (83) 

fe=i fcc * Z fc =i 

and 

1 oo 

(a S;S |H|a S;S ) = ft^|a s ^| 2 + -Tiuj r E kc l\fk(s, k)\ 2 - (84) 



2 

fe=i 



The two averages will differ only by the value of the vacuum contribution if Cfe = 1/ \fk which leads back to the above 
mentioned choice of c^. With this choice and taking the general combination of coherent states 

|*) =]>> s , b Kk) (85) 



wc find 



CO 

= e ^i<M 2 Kd 2 + 9 E ^i<m 2 E fc i/ fc ( s > «)i 2 ( 86 ) 



2 

S.K fc — 1 



and 



(¥|H|¥) = E ^«|^| 2 |a^| 2 + J E ^l^.-l 2 - ( 87 ) 
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B. Vacuum 

Similarly to the one-oscillator case the traditional notion of a vacuum state is replaced in our formalism by a vacuum 
subspace consisting of all the vectors of the form 

l*) = E Vl/ il,ol^«A,0) 

® E *i 2 !a 2 ,«x 1 ,«x 2 ,0,ol S l'«A 1 ,0)|s 2 ,K A2 ,0) 

® E *S, 2 ,,3^ 1 ^ 2 ^ 3 AO,ol S l'«A 1 ,0)|s 2 ,K A2 ,0)| S3 ,KA3,0) 

© . . . (88) 

It seems that there is no reason for introducing the standard "vacuum state" understood as the cyclic vector of the 
GNS construction. 

In the discussion of various vacuum phenomena (e.g. spontaneous emission) we will assume for simplicity that all 
the oscillators are "embedded" in identical vacua i.e. the multi-oscillator vacuum is of the form 

®Vpi\4>)\4>) 

®^\<t>)\<t>)\<t>) 

© . . . (89) 
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The average energy of the free-field vacuum state is therefore 

oc 

H= =J2 n Pn(<f>\H\<l>) =nH (90) 

n=l 

where n and H are, respectively, the average number of oscillators and the average energy of a single oscillator. Again 
no problem with infinite vacuum energy is found. Obviously, one can contemplate also other vacua, say, in entangled 
or mixed states. 



C. Multi-oscillator vs multi-photon states 

The coherent states we have introduced at the one-oscillator level involve superpositions of different excited states. 
We know that in the traditional approach the transition between two such states, say, 

|«,/6,2)->|s,*,0) (91) 

is interpreted as an absorbtion (by some system) of two photons. In the new formulation the problem is more 
complicated since the "2-photon" absorbtion may be represented also by 

\s, Ft, l)\s, Ft, 1) -> \s, Ft, 0)|s, Ft, 0). (92) 

The two types of transitions do not represent the same process and the two final states are physically distinguishable 
since their energies are different. Indeed, 

H\s, Ft, 0) = H\s, Ft, 0) = ^huj R \s, Ft, 0) (93) 

whereas 

n\s,Ft,0)\s,Ft,0) = (H®l + l®H)\s,Ft,0)\s,K,0) 

= huj H \s,K,0)\s,Ft,0). (94) 

The notion of a 2-photon state becomes therefore somewhat ambiguous. To make it more precise one has to formulate 
a photodetection theory within the new framework. In what follows we shall try to avoid the use of the word "photon" 
and will talk about "light quanta" and "multi-oscillator" (or n-oscillator) and "higher-excited" (or n-th excited) states 
of light. 

The 2-oscillator states 

|*±>= XJ 4 2 1 ) ,« 2 ,„(l+^i^)|-^2,n)±|-,K 1 ,n)|+,K 2 ,n)), (95) 

satisfying 

* ( - 2) - =±* ( - 2) - (96) 

are (for any n > 0) perfectly justified generalizations of the standard 2-photon maximally entangled state. We shall 
later see that although such "higher excited photons" (i.e. n > 1) are in principle possible, they are not produced in 
a two-photon spontaneous emission (at least up to second-order perturbative effects). The technical reason for this is 
the same as in the ordinary formalism and is related to properties of the annihilation operator a. 



VII. SPONTANEOUS EMISSION OF A "SINGLE PHOTON 



In this section we shall again consider the spontaneous emission of light within the two-level-atom approximation. 
The example illustrates some pecularities of the multi-oscillator formulation. 

Denote by Hf the multi-oscillator Hamiltonian of the free field we have discussed in the previous two sections. The 
dipolc and RWA Hamiltonian of the 2-level atom interacting with quantized electromagnetic field is now 
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H= ^ttij a 3 +H F + huodY^ {gs,K X a St u x o+ + 5*^0* ^cr_). (97) 

Similarly to the one-oscillator case one has 

e iH F t/h aa ^ e -iH F t/h = e -™*t asR (98) 

and therefore the interaction-picture Hamiltonian is 

Hi = hw dJ2 (g s ^e i ^-^ t a s ^a + + g^e^^al^) . (99) 

S,K X 

The first pecularity we encounter is the fact that the Hamiltonian is block diagonal with respect to © and therefore does 
not have nonvanishing matrix elements between spaces corresponding to different numbers of oscillators. As a result 
the interaction cannot change the number of oscillators, a property of crucial importance for statistical properties of 
light as we shall see in the context of the Planck blackbody radiation law. 
Consider the initial state 

l*(0)> = \/Pil+>^<Ms,«,0) 

S,K 

l,Kl -^2 

Si ,Ki S 2 ,K2 

®VP3\+) E E ^ 1 ^ 1 S2 ,K 2 ^ S3 ^3l S l^l. O >l S 2,K2,O)|s 3 ,K 3 ,O) (100) 

Sl,Kl S 2 .K2 

The first-order perturbation theory yields (it is instructive to keep again the constants c„ arbitrary) 

-i(u> -u>z)t _ 1 

|*(t)> = |*(0)) + ujodc^Y 9U4>sM-)\s, *, 1) 

S,K 

( g-iCwo-waJt _ j e -i(wo-u>B 2 )t _ j \ 

©Wodc 2 Vp^|-> 2^ ff; i ,a 1 ^«i,Bil*i.'«i.l>l0>+ 2^ 5 s * 2 ^ 2 ^ 2 ,K 2 |0)|s 2 ,K 2 ,l) 

\ Sl ,Kl S 2 ,K 2 / 

©... (101) 

As we can see, the "single-photon" emission can be realized in an infinite number of different ways. In the 1-oscillator 
subspace the oscillator simply gets excited to the 1-st excited state. The probability amplitude for this process is 
proportional to the probability amplitude that the field is found in a 1-oscillator state. In the 2-oscillator subspace 
there are two possibilities: Either the first oscillator gets excited and the second one remains in the ground state, or 
the other way around. The probability amplitude for this process is proportional to to the probability amplitude that 
the field is found in a 2-oscillator state. And so on. 

Repeating the argument given for a single-oscillator description, assuming the isotropy and polarization- 
independence of the vacuum mode density, we arrive at the spontaneous emission rate of the form 

00 

P = FMPoid E nc «P« ( 102 ) 

71=1 

where F(un) = \4> s m\ 2 - As we can see the choice c n = 1/v 7 "- plays again a special role since then 

P = F(w )P old , (103) 
that is, the result is the same as in the single-oscillator description. 

VIII. SPONTANEOUS EMISSION OF "TWO PHOTONS" 

In what follows we will start with the Hamiltonian H = H + V, where 
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H A + H F (104) 

A(x) ■ p 

m 

i?A is the full (i.e. infinite- level) Hamiltonian describing an atom and Hi? is the free-field Hamiltonian dicussed in 
Sec. V and obtained by the multi-oscillator extension of the one-oscillator Hamiltonian introduced in Sec. III. To 
simplify notation we shall denote the sum and the integral over, respectively the discrete and the continuous parts of 
the spectrum of Ha by the sum ^2 C - We are not making the rotating wave approximation. In the dipole approximation 
we set x — 0. We shall also keep the constants c n arbitrary. 



#0 = 

V = 



A. Two different light-quanta in 2-oscillator subspace 

In this subsection we will use the second-order perturbation theory to compute the amplitude 

(6|(s 1 ,K 1 ,l|(s 2) K 2 ,l|Z7(t / ,t i )|a)|*) (106) 

where the states |sfe,Kfe,l), k — 1,2, are orthogonal, U(tf,ti) is the evolution operator mapping the initial state at 
time U into the final state at time f/, |^) is a vacuum state (p9|), and |a), \b) are two bound states of the atomic 
Hamiltonian Ha- 

The fact that the interaction term does not change the number of oscillators reduces the above amplitude to its 
2-oscillator counterpart 

p^(b\(s 1 ,K 1 ,l\(s2,K 2: l\U(t f ,t l )\a)\^) (107) 



Using the standard perturbative techniques we obtain the second-order approximation H (see Appendix XI C) 

(&| (S 1 , «1 , 1 1 <S 2 , K 2 , 1 1 ^ (2) (*/, *i ) k> I 1 0) 



2 2nie 2 I h I h \- (e7 2 ,g 2 ' Pbc) (e* ^ • Pea) (T) 



^K 2 1 



2 27T Je 2 / ft / h \- (e a * , gl ■ ^c) (g* ^ ; Pea) (T) 

V 2^,1/ V 2w Sl K ^ K,«i,o,« 2 ,o - ^c,«i,o,« 2 ,i + «0+ 



where <5( T )(£ - £') = [tt(J5 - i?')] _1 sin ((E - E')T/2h),p bc = (b\p\c), and p ca = (c\p\a). 
The energies occuring in the above expression are (ground-state energies are not removed!) 

■Ea,*i,0,*2,0 = ^ + T^^i + ^ w «2 ( 109 ) 

I = -Eft + (1 + 7;)^ + (1 + 2)^2 ( 110 ) 

/•'<,• .I,;,,. = Be + (1 + 2)^1 + \ huj Z2 (HI) 

/•'.,.- 1 = B c + 2^ gl + (1 + i)7k^ 2 (112) 

The net result is the following 

(b\(s U K U l\(s2,K2,l\U(t f ,ti)\a)\y) 



(e * p*b ) (e * p ) 



.2-Kie 2 



(e * pi, ) (e* p ) 

Let us note that the amplitude is symmetric with respect to permutation of states of the two oscillators: 
(b\(sx,Kul\(s 2 ,K2,l\U(tf,ti)\a)\y) = (6|<a a , >?a, l|<*x, ^i, **)l a >l*> 
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B. Two different light-quanta in 3-oscillator subspace 



Consider the amplitudes 

(6|<si, ^i, 1]<*2, i^, lK^ll/Ct/,^)^)!*) (113) 

(6| <S X , , 1 1 <01 <* 2s /^!, 1 1 *i) I *»> I ( 114 ) 

(b\{<t>\{si,Hi,l\{s 2 ,&A\U(tf,U)\a)\*) (115) 

In the framework we propose it is necessary to include the contributions of this type arising from all the possible 
numbers of oscillators. 

It is again sufficient to restrict the analysis to 



^{b\{s 1 ,K 1 ,l\{s 2 ,K 2 AMU(tf,ti)\a)\4>)\4>)\4>) (116) 
and similarly with the other two amplitudes. In second-order perturbation theory (see Appendix XI D| ) 



c^ 2 <6K Sl ,K 1 ,lKs2,K 2! lK0|C/ (2) (t/,ti)|a)|^)|^)|^) =% 2 (b\(s 1 ,K 1 ^\(s 2 ,K 2 A\uW(t f ,t i )\a)\(i>)\<t>)- (117) 

The result is therefore essentially identical to the one obtained for the 2-oscillator subspace. 

A closer look at the derivation of the 3-oscillator contribution shows that (i) exactly the same will happen for any 
number of oscillators and (ii) the second-order amplitude describes an emission of at most two quanta. 

C. Comparison with the standard formalism 

It is instructive to compare the result we have obtained with the second-order calculation performed by means of 
the ordinary quantum optics formalism. Let 

\s 1 ,ki;s 2 ,k 2 ) = \s 2 ,k 2 ;s 1 ,k 1 ) = aJ i)Si aJ 2>g2 |0) (118) 

be the two-photon state of the standard formalism, a 1 - the standard creation operator and |0) the vacuum state. 
Then 

{b\{si,Ki;s 2 ,K 2 \U^(t f ,U)\a)\0) 



2nie 2 v ^ / ft / ft _ ^,_ M , _ > £ (T) (£ a - E b - huj Kl - ftuj K2 ) 

■ m? ^\ 2u, Xl V\ 2u Ka V (e '^'^ e '>*' Pea > E a -E c -ftu> K2 +iO + 



2-Kie 2 \ / ft I ft _ _ v S^^Eg - E b - ftuj Kl - ftui K2 ) 

m 2 ^V^^V^^ E a -E c -hv K1 +iQ+ ■ 

This is precisely the same expression that occurs in the modified amplitudes. 

D. Probability of spontaneous emission of two quanta 

In the subspace corresponding to n oscillators the "two-photon" emission can take place in 

n ^ _ n(n - 1) 
2 



(119) 



different ways. Taking into account probability amplitudes associated with all the n-oscillator subspaces, n > 1, and 
the fact that the two quanta can be emitted in two different orders, we obtain 

p(s 1 ,k 1 ,s 2 ,k 2 ) = 5 c t,Pn\<l>s 1 ,it 1 \ 2 \<l>s 2 ,ii 2 \ 2 p(si ) Kl, S 2 , K 2 ) oM (120) 

n=2 

where p(sx, Hi, s 2 , K 2 ) \d is the standard result obtained by means of ordinary quantum optics. Taking, as before, 
c n = 1/v 7 " wc nn( l 
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oo 

p(si,Ki 1 S 2 ,K 2 ) = (l - X! ~ ) \ ( I ) S 1 ,K 1 \ 2 \4>S 2 ,K 2 \ 2 P(S1 1 Kx, S 2 , K 2 )old (121) 

n=l U 

Under such assumptions the angular distribution of the two-photon emission is the same as in the standard theory. The 
probability of the 2-photon spontaneous emission is thus a product of four terms. It may be difficult to distinguish 
between i ;l (w^ 1 )F(co'^ 2 )(l — (-)) and analogous factors arising in real experiments from detector inefficiency. The 
above result may have therefore nontrivial implications for the problem of testing quantum mechanics versus local 
hidden-variables theories and is very closely related to the so-called detector inefficiency loophole in Bell's theorem 
(cf. flO|jrT|| ) . The reason is that the presence of F(w^ 1 )F(wg 2 )(l — (— )) will necessarily lower two-photon coincidence 
rates, whereas it is known that in order to violate the Bell inequality the rates must exceed certain thresholds. The 
problem is worth further studies. 



IX. BLACKBODY RADIATION 



One of the possible tests of the new formalism is the problem of blackbody radiation. Planck's famous formula |T^] 

where N w is the average number of excitations of an oscillator in inverse temperature (3, is one of the first great 
sucesses of quantum radiation theory and marks the beginning of quantum mechanics. Contemporary measurements 
of q(uj) |l^,Q performed by means of COBE (Cosmic Background Explorer) are in a very good agreement with 
the Planck law. The data have been carefully analyzed in the context of nonextensive statistics (l^Jl^l in search of 
possible deviations from extensivity. The result that comes out systematically is \q — 1| < 10~ 4 where q is the Tsallis 
parameter. The case q — 1 corresponds to the exact Planck formula. If there are any corrections whatever, they must 
be quite small. 

The standard derivation of the formula consists basically of two steps. First, one counts the number of different 
wave vectors k such that c\k\ € [lo,lo + Aw], Second, one associates with each such a vector an oscillator and 
counts the average number of its excitations assuming the Boltzmann-Gibbs probability distribution at temperature 
T and chemical potential /i = 0. The latter assumption is justified by the fact that the number of excitations of the 
electromagnetic field is not conserved in atom- light interactions. 

In the new model the situation is slightly different since there exists an additional conserved quantum number: 
The number of oscillators . As we have seen in previous calculations the Hamiltonian is block-diagonal with respect 
to © but changes the number of excitations in each ./V-oscillator subspace of the direct sum. The state vectors at 
the multi-oscillator level are symmetric with respect to permutations of the oscillators and therefore the oscillators 
themselves have to be regarded as bosons whose number is conserved and their chemical potential is ^ ^ 0. However, 
their excitations should be regarded as bosons with vanishing chemical potential. 

The energy eigenvalues 



>(n + \) (123) 



corresponding to the oscillator whose frequency is w are parametrized by two natural numbers: m (the number of 
oscillators) and n (the number of excitations). Assuming the standard Boltzmann-Gibbs statistics we obtain the 
probabilities 

p mn = £-i e -/3[mM«+l)-«v] (124) 

where 

Z=Y e /MP+iWa)_f . (125) 

m— 1 



The Lambert series |17 



1^ «™ — 



(126) 
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is convergent for any x if J2 m =i Qm ^ s convergent. Otherwise (126) converges for exactly those x for which the power 
series YZ=i a ™x m does. In (|l2f) a m = e /9™0*+ftu>/2) and YZ=i a m < oo if fx + fkj/2 < 0. If p, + huj/2 > we still 
have convergence of (125) as long as X)m=i er^ m ^ n ^~^ < oo. The upper limit imposed on fi by the finiteness of Z 



is therefore ^Tiuj — /x > 0. In what follows we assume that fi is w-independent and therefore fi < 0. 
The appropriate average number of excitations is 



mne-^ mRuJ{n+ ^ ) - m ^ (127) 

m— 1 n— 

and the Planck formula is replaced by 

fticwM = -^rw 3 n w . (128) 

It is easy to show that £> n ew(w) tends to the Planck distribution with /i — > — oo. To see this consider a more general 
series 

oo oo 

Z- 1 ^ mnq m e-' 3mhuj( - n+ ^ (129) 

m— 1 n— 



-/3[m^tj(n+4)— m/x] 



where Z is the normalization factor and X)m=i Qm < oo- If Si = 1 and g m = for m > 1 then ( 129 ) is just the exact 
Planckian formula. Factoring out e - ' 9 '' 1 ' in both the numerator and the denominator of n w we obtain qi = 1 and 
q m = e - ' 3 ''^™ -1 ) for m > 1. For — > oo all q m , for to > 1, vanish and the limiting distribution is Planckian. 

This proves that an experimental agreement with the ordinary Planck's q(uj) cannot rule out our modification 
but can, at most, set a lower bound on an admissible value of However, assuming that /i has some finite 
and fixed value it should be in principle measurable. The plots show that the modifications become visible around 
(i sa —3ksT. Assuming that the chemical potential is temperature independent, say // = —ksTo, we obtain a kind 
of critical temperature T cr i t i ca i ~ To/3 above which the ratio p/iksT) is small enough to make the modifications of 
the distribution observable. For T < Tcriticai the distribution should be given by the Planck law; for T > T cr iticai the 
distribution should approach the /i = distribution, i.e. this would be a Planck-type curve but with the maximum 
lowered and shifted towards higher energies. 

Fig. 1 shows the plots of p new (w) for p — (lower dotted), fi — — 0.8fcsT (upper dotted), and /i = — lOfc^T (solid). 
The thick dashed curve is the Planck distribution. The curve obtained for (i — —lOkgT is indistinguishable from the 
Planck distribution. The plot does not change if one takes fi < — lOfc^T and differences are not visible even if one 
plots the distributions in logarithmic scales (not shown here) . This is a numerical proof that the distribution we have 
obtained on the basis of the modified quantization tends very quickly to the Planck one as (i — > — oo. It is instructive 
to compare the modification we have predicted with those arising from nonextensive statistics. The two thin dashed 
lines represent Tsallis distributions resulting from nonextensive formalism for q = 0.95 (lower) and q = 1.05 (upper). 
The modifications we have derived are therefore qualitatively different from those resulting from Tsallis statistics. 



\ 



\ 
\ 
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FIG. 1. Q ncv ,(u>) for jj, = (lower dotted), \i = — 0.8/csT (upper dotted), and /i = — WksT (solid). The energy range 
is O.OlfcsT < Thti < lOfcsT. The thick dashed curve is the Planck distribution. The curve obtained for fj, = — lOfcsT is 
indistinguishable from the Planck distribution. The two thin dashed lines represent Tsallis distributions resulting from the 
Tsallis formalism for q = 0.95 (lower) and q = 1.05 (upper). Since g ncv < g at least in the neighborhood of the maximum, 
the new distribution has to be compared with q < 1 statistics. The curves are qualitatively different. In particular, all q < 1 
distributions require an energy cut-off which does not occur for £» n cw 




2 4 6 8 10 

FIG. 3. Contour plot of £> n ew(w) for — lOfcsT < fj, < 0. The fast convergence to Planckian q(uj) (as n — * —oo) is clearly 
seen. 
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X. CONCLUSIONS 



"A theory that is as spectacularly successful as quantum electrodynamics has to be more or less correct, although 
we may not be formulating it in just the right way" . The above quotation from Weinberg could serve as a motto 
opening our paper. The main idea we have tried to advocate was that the standard canonical quantization procedure 
is, in certain sense, too classical to be good. 

The reasons for such a choice of quantization could be both historical and sociological and may be rooted in the fact 
that the idea of quantizing the field was formulated before the real development of modern quantum mechanics. In 
oscillations of a simple pendulum it may be justified to treat ui as an external parameter defining the system (via, say, 
the length of the pendulum). But oscillations of the electromagnetic field do not seem to have such a "mechanical" 
origin and it is more natural to think of the spectrum of frequencies as eigenvalues of some Hamiltonian. That is 
exactly what happens with other quantum wave equations. 

We have defined the quantum electromagnetic field as an oscillator that can exist in a superposition of different 
frequencies (or, rather, wave vectors). This should not be confused with the classical superpositions of frequencies 
created by, say, a guitar string. The superpositions we have in mind dissappear at the classical level. 

Once one accepts this viewpoint it becomes clear how to quantize the field at the level of a single oscillator. We 
do not need many oscillators to perform the field quantization. But there is no reason to believe that all the possible 
fields can be described by the same single oscillator. And even more: We know that the structure of the one-oscillator 
Hilbert space is not rich enough to describe multi-particle entangled states and there is no doubt that such states are 
physical. The next step, performed already after the quantization, is to consider fields consisting of 1, 2, 3 and more 
oscillators, and even existing in superpositions of different numbers of them. The resulting structure is analogous to 
the Fock space but so the procedure can be (although somewhat misleadingly) referred to as "second quantization" . 
What is essential we do not need the vacuum state understood as the cyclic vector of the GNS construction. In the 
new framework such an object seems rather artificial. 

On the other hand, there exist vacuum states. These are all the states describing ground states of the oscillators. 
They correspond to concrete finite average values of energy. A general vacuum state is therefore a superposition of 
different eigenstates of a free Hamiltonian and is not, in itself, an eigenstate of the Hamiltonian. 

No further assumptions are made. The system is described by laws of ordinary quantum mechanics so that to 
compute concrete problems we can use standard methods. Perturbation theory leads to structures we know from the 
standard Feynmann diagrams. The blackbody radiation is calculated by means of the standard Boltzmann-Gibbs 
statistics. 

Let us close these remarks with another quotation: "Present quantum electrodynamics contains many very impor- 
tant 'elements of truth', but also some clear 'elements of nonsense'. Because of the divergences and ambiguities, there 
is general agreement that a rather deep modification of the theory is needed, but in some forty years of theoretical 
work, nobody has seen how to disentengle the truth from the nonsense. In such a situation, one needs more exper- 
imental evidence, but during that same forty years we have found no clues from the laboratory as to what specific 
features of QED might be modified. Even worse, in the absence of any alternative theory whose predictions differ 
from those of QED in known ways, we have no criterion telling us which experiments would be relevant ones to try. 

It seems useful, then, to examine the various disturbing features of QED, which give rise to mathematical or concep- 
tual difficulties, to ask whether present empirical evidence demands their presence, and to explore the consequences 
of the modified (although perhaps rather crude and incomplete) theories in which these features are removed. Any 
difference between the predictions of QED and some alternative theory, corresponds to an experiment which might 
distinguish between them; if it appears untried but feasible, then we have the opportunity to subject QED to a new 
test in which we know just what to look for, and which we would be very unlikely to think of without the alternative 
theory. For this purpose, the alternative theory need not be worked out as completely as QED; it is sufficient if we 
know in what way their predictions will differ in the area of interest. Nor does the alternative theory need to be free 
of defects in all other respects; for if experiment should show that it contains just a single 'element of truth' that is 
not in QED, then the alternative theory will have served its purpose; we would have the long-missing clue showing in 
what way QED must be modified, and electrodynamics (and, I suspect, much more of theoretical physics along with 
it) could get moving again " |3j. 
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XI. APPENDIX: TECHNICAL DIFFERENCES AND SIMILARITIES WITH RESPECT TO THE 

STANDARD FORMALISM 



A. Proof of Eq. (11) 

This calculation is elementary but very important, so we give it explicitly: 

a Uk (t) =e m l n a^ m l h (130) 
= e ifi®(ata+|l)t| a;fe ^ Wfe | 8ae -iQ®(ata+il)t 

= |w fc )(w fc | (8 e *-*(« t «+4i)*ae-* , '*( ato+ * 1 )* (132) 

= \w k ) (w fc | ® e^a = e-*"**^ (0) (133) 



B. Energy-momentum operators for free fields: 1-oscillator formulas 

To see how (^) and (0) arise let us first note that 

a Sl ,K l a S2 .R 2 = (134) 

a\ g a\ g = (136) 

„„J /„t \2 



unless si — s 2 and «i = k 2 [cf. Eqs. ( |27|) and (|2§)]< The terms involving (a Sj s) 2 and (a' s -) 2 disappear due to 
e s ,re • e s s = and its complex conjugate. As a result 



- (f (i, x) ■ k{t, x) + 3(jt, x) ■ 3{t, £)) = e lP - x ' n ^2^\s, k x ) (a, «a| ® + 0*0) e^"*/* 

S,K A 

= XI ^tI s ' k a)(s,k a | $ (aa^ + a^a) = H/V. 

To find the relation between the Pointing vector and linear momentum we first have to show that 

E x B = -B x E. (137) 

The relevant formula is 

[E a ,Bp] = ^2 itw\-z[8a0 - (nit x ) a (riK X )fij\s,Kx){s, K\\ <8> 1. (138) 



2 

The remaining calculations are similar to those for H. 

C. Derivation of the "2-photon" amplitude: 2 oscillators 

We employ the standard second-order time dependent perturbation theory and notation from 
c^ 2 <6K«i,/5i,l|(« 2 ,/52,l|^ (2) (t/,ti)|a)|^)|0) 

27rie 2 ^ ^ <5' T '(-^a-fti,0,K2,0 — ^6,«i,l,re 2 ,l) 
= ~^" 2-, ^r.fc^r'J 1 ^ ; 



'' k ' E j £t + i0 + 

c, Si, S2,-Ki,if2,rai, "2 r,k,r',k' L ' 11 " 



(6|(si, Ki, l|(s 2 , k 2 , 1| X J ^(a s ,se s ,j • p + a^e*K ' p) l c ) |5"i, n-i) \S 2 ,K 2 , n 2 ) 
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x (c\{S 1 ,K u n 1 \{S2,K2,n 2 \ J2 \ 2^v( as '^' gs '^' ' ^ + a ^'^'<* ' p)\ a )\ r ^,0)\r' ,k' ,0) 
^irie 2 A ^ ^^ T H-^o,«i,o,« 2 ,o - ^6,rei,i,re 2 ,i) 

m 2 Z-^ 



r'k' 



'* k ' E r „7, „ — E z £t + iO + 

c,Si,S2-Ki,K2,n-i,ri2 r,k,r' ,k' ' ' ' ' 



(si,Ki, 1|(s 2 ,k 2 , 2cJ~X^ ^ 1 + 1 ® a s,a)esM-Pbc + (a\ R ® 1 + 1® a\ R )e* R ■ p bc 



x \Si,K 1 ,n 1 )\S 2 ,K 2 ,n 2 )(S 1 ,Ki, ni\(S 2 , K 2 ,n 2 



^ V 2uj i -,v i^ as ''' t ' 1 + 1 <g)as '.«')'^'^' '^a + «) 1 + 1 Oal',K')e*',K' • Pea) |r, fc, 0)|r', fc', 0). 



x 

S' ,K 



The block-diagonal property of the interaction Hamiltonian has been used twice. The remaining calculations arc 
standard. It is remarkable that although the result we obtain is essentially the same as in the standard formalism, 
the technical reasons for this are completely different. 

D. Derivation of the "2-photon" amplitude: 3 oscillators 

Here we sketch the proof of the 3-oscillator amplitude. In second-order perturbation theory 

C 3 2 <6j< Sl ,« 1 , i|< S2 ,« 2 ,iK0l^ (2) (*/>*i)l«>l0>!0)!^> 

c,Si)S2,Ki,K2,K3,Tiii'ri2,'n3 ro,ko,r,k,r f ,k' ,r" .k" 



x (6|(si,Ki,l|(s2,K2,l|(ro,fco,0|^ J -^-y(a s #e s # -p + a\ -e*« ■ P^\c)\Si, Ki,n 1 )\S 2 , K 2 ,n 2 )\S 3 , K 3 ,n 3 ) 



S.K 



x (c|(S 1 ,^ 1 ,n 1 KS 2 ,^ 2 ,n 2 |(53,^3,n 3 |^W^^(a 5 ,^e s ,^ • p + a\, R ,e* s , ~ K , ■ p) |o)|r, k, 0}|r', k', 0)|r", k", 0} 

^^ T '(-^a,fe,0,fc',0,fc",0 _ ^6,Ki,l,« 2 ,l,feo,o) 



^o,fc,0,fc',0,fe",0 ^c,K 1 ,n 1 ,K 2 ,n2,K 3 ,n 3 + *0+ 

c,Si,S2,-f^i,^2)^3,7ii,Ti2)^3 r ,k ,r,k : r' ,k' ,r" .k" 



x (6|(si,Ki, l|(s 2 ,K 2 ,l|(r ,fc ,0| ^ ^ ( a ^,«^,« ' P + a Igg s *K ' -P 



x | C )|5 1 ,^ 1 ,n 1 )|5 2 ,^ 2 ,n 2 )|,S3,^3,n3)(c|(5 1 ,K 1 ,n 1 |(5 2 ,if 2 ,n 2 |(53 ) if3,n3| 



K 2^ \/ ^7 ( a ^«' ® 1 ® 1 + 1 ® 4,* ® 1 + 1 ® 1 ® a\,^e* s ,^ ■ p\a)\r, k, 0)|r', fc', 0) |r", fc", 0) 

s' ,K.' V 

J \ O,fe,0,fe',0,fe",0 6,Kl,l,K2,l,fco,0/ 



^o,fe,0,fe',0,fe",0 c,Jfi,ni,K 2 ,n2,^3,n3 + 



The remaining part of the proof is standard. In the course of the computation one recognizes the elements known 
from standard Feynman diagrams, in particular the self-energy corrections due to emission and reabsorbtion of virtual 
photons. A general property of the perturbation series we find is its better convergence due to the presence of the 
vacuum amplitudes <f) s g. 
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